WIS99/27/Aug-DPP 



|hep-th/ 0005075 



On the Bosonic Spectrum of QCD(1+1) with 

SU(N) Currents 



Uwe Trittmann 

Department of Physics 
Ohio State University 
Columbus, OH 43210, USA 



Abstract 



In this note we calculate the spectrum of two-dimensional QCD. We formulate the 
theory with SU(N C ) currents rather than with fermionic operators. We construct the 
Hamiltonian matrix in DLCQ formulation as a function of the harmonic resolution 
K and the numbers of flavors Nf and colors N c . The resulting numerical eigenvalue 
spectrum is free from trivial multi-particle states which obscured previous results. The 
well-known 't Hooft and large Nf spectra are reproduced. In the case of adjoint fermions 
we present some new results. 



1 Introduction 



On the way to understand the physics of strong interactions, two-dimensional QCD 
has remained an interesting model ever since it has been studied by 't Hooft [![. It 
has, however, the disadvantage of lacking dynamical (i.e. transverse) gluonic degrees 
of freedom. To construct a model closer to four dimensional QCD, adjoint rather 
than fundamental fermions were built into the theory Also matter with a large 

number of flavors has been considered [|J || . Surprisingly, these theories are related by 
a universality ||: the massive spectrum and interactions of two-dimensional Yang-Mills 
theories coupled to massless fermions in arbitrary representations depend only on the 
gauge group and the level of the afhne Lie algebra. The associated parameter space of 
two-dimensional Yang-Mills theories with massless fermions is depicted in Fig. [I] 0. 

Of particular interest in two-dimensional QCD is the transition between confinement 
and screening. The 't Hooft model is known to consist of stable mesons and has no 
confinement-screening transition. The model with adjoint fermions exhibits deconfine- 
ment at zero fermion mass m, and the string tension vanishes like o oc m at this point 
HI II GU • Interestingly, this theory is supersymmetric at the fermion mass m? = g 2 N c , 
and asymptotically supersymmetric otherwise. It has been studied extensively in the 
literature 0, |], [TT|, [12], |13| . We will focus on the massless case of this theory, where 



we can make use of the universality and consequently we will formulate the theory 
as a current algebra problem. The Hilbert space of the problem then splits up into 
current sectors. This allows for a convenient reduction of the numerical effort when 
solving for the spectrum of the theory, since all single-particle states lie in the current 
block of the identity (bosons) or in the adjoint current block (fermions) ||. Because 
of its simpler algebraic properties, we will limit ourselves in the present work to the 
bosonic spectrum. We will use the additional parameter in the problem, A = Nf/N c , 
as a tool for interpreting the spectrum, i.e. we will consider theories 'on the arc' of 
Fig. |l|, where both N c and Nf are large. The main emphasis will be on the adjoint 
case (N c =Nf). The goal is to identify the single-particle states of the theories and to 
be able to connect to results anticipated from the asymptotic theory, e.g. the expected 
multi-Regge trajectory structure in the adjoint case ||. 

A very convenient way to formulate two-dimensional QCD is to quantize the system 
on the light-cone |]J. Usually the light-cone gauge, A + =0, is used. This approach gives 
a very complete picture and includes non-perturbative effects, if an appropriate regu- 
larization for the quark self-energy is used ||14|| . Discretized Light Cone Quantization 
(DLCQ) [|15| is then a method especially suited to solve numerically for the spectrum 
of low dimensional theories. The momenta are discretized by imposing boundary con- 
ditions on the fields. The typical DLCQ program is to construct a finite dimensional 
Hamiltonian matrix characterized by the harmonic resolution K. The spectrum is ob- 
tained by diagonalizing this matrix numerically for larger and larger K, and eventually 
extrapolating to the continuum limit, K —>■ oo. 



The latest large N c analysis of adjoint QCD 2 [TT[ revealed several single-particle 
states and, most interestingly, a continuum of states at precisely four times the mass 
(squared) of the lowest fermion state. Despite this remarkable structure, the insight 
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Figure 1: The parameter space of two-dimensional Yang-Mills theories with massless 
fermions. 



gained from this finding is rather small. It is clearly not a signature for screening versus 
confinement. Both in a confining and a screening theory, one would expect a continuum 
of states at 4m 2 of a single-particle state, although from very different mechanisms. 
More information about the spectrum is thus needed to understand things like the 
deconfinement mechanism or the statistics of single-particle states. 

The paper is organized as follows. In the following section we review two-dimensional 
QCD with massless fermions. Sec. |3| deals with analytic considerations. Namely, we 
show that it is possible to calculate two eigenvalues of the theory from the outset. In 
Sec. |4.1| we construct the Hamiltonian matrix for any harmonic resolution K and arbi- 
trary {i.e. in principle finite) Nf and N c . We will, however, in praxi only use the large 
N limit of this result. Sec. |] details the numerical results. We will focus on the adjoint 
case, but also consider the 't Hooft limit and the large Nf limit to present a coherent 
picture of two-dimensional QCD. A discussion of the results follows. 



2 QCD in two dimensions 

We want to compute the massive spectrum of SU(N C ) Yang-Mills gauge fields coupled 
to massless adjoint fermions in two dimensions. Due to the universality in these gauge 
theories 0, this is equivalent to solving for the massive spectrum of Nf=N c flavors 
of massless fundamental Dirac fermions coupled to the gauge fields. The latter case 
turns out to be more general, in the sense that we have an additional continuous 
parameter, X:=Nf/N c , in the theory. The case of adjoint gauged fermions is recovered 
when setting A to unity. Without loss of generality we consider only the Lagrangian of 
the fundamental theory 

i N f _ 

£ = Tr[-—F, U F^ + i £ * a% D^ a } (1) 

where ^5f a = 2 _1//4 (^ a ), with ip a and Xa carrying color indices, which we suppressed. 
The field strength is F^ v = d^A u — d u A^ + £[.A M ,A„], and the covariant derivative is 



defined as = + i[A^, ■}. We work in the light-cone gauge A + = and use the 
convenient Dirac basis 7 = a\, 7 1 = — ioi. The Lagrangian then becomes 



C = Tr 



—r{d-A + f + # t «9+^ + ix ] d-x ~ A+J 



2g 2 



(2) 



with the current = ijj\ a ip a j. We can integrate out the (non-dynamical) component 
A + of the gauge field and obtain 



£ = Tr 



(3) 



It is obvious that the left-movers x decouple. Noting the simple expression of the 
interaction in terms of the currents, it is natural to bosonize the theory [|16| . We follow 
Ref. || to derive the momentum operators. The bosonized action of colored flavored 
fermions is 

So = S^ z /(g) + S^ w (h) + l -J d 2 xd^dy } (4) 



where g G SU(N C ), h G SU(Nf), exp{i^J N j^ 4>} G Ub(1), and the Wess-Zumino- 
Witten action is 



qWZW 



(g) = £-[ d 2 xTr(d,gd»g- 1 ) + -^ 

87T JS lZ7T JB 



(Pye^Trig^d^g-'d^g-^g), (5) 



where B is the solid sphere whose boundary £ represents space-time. The action of 
the theory is then 



(6) 



with the currents J = 7^gd_g 1 generating a level Nf affine SU(N C ) Kac-Moody 
algebra. The associated energy-stress tensor T^ u yields the momentum operators 



P A 
P- 



7T 



N c + N 
,2 



(x : 



(7) 



g l roo I 
T = — / dx . J ; L j{x ) — 2" Jji{x ) . . 
Z J — 00 u_ 



To obtain the mass eigenvalues M n we have to solve the eigenvalue problem 

2P+P-\?P) = M^), (9) 

which is equivalent to diagonalizing the operator P~, since P + is already diagonal. 
To discretize the system, we impose periodic boundary conditions of length 2L on the 
currents, Jij(x~ — L) = Jij(x~ + L), and expand them into a discrete series of modes 



Jij («£ 



1 



A' 



V2L 



n e 



(10) 



n=-K 
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with 

Jy(n)|0) =0 Vn>0. (11) 

The cutoff K = P + L/2n controls the coarseness of the momentum-space discretization. 
The continuum limit is obtained by sending K to infinity. The modes of the currents 
obey the algebra 

[Jij(n), Ji'f(n')] = nNf (SiySvj - j-j-bijSi'j^j $l n ,+6ij>Ji>j(n+n')-6iijJiji(n+n'). (12) 
The momentum generators 

P+ = (t)a£T1v> t Jii( - n)Unl (13) 

p ~ = f; I : ^M-n)JiiM, (14) 

are finite-dimensional matrices on the Hilbert space constructed by acting with the 
current operators on the vacuum defined in Eq. ([11]). For convenience we introduced 
the scaled coupling g 2 = ^— . Note that the box length L drops out of the eigenvalue 
problem, Eq. ([5]). At large N c we expect the Fock basis to consist of single-trace states 

±-Tr[J(-n 1 )J(-n 2 )---J(-n s )]\0). 

c 

We note that a cyclic permutation of the currents will reproduce these states only up 
to states with a lower number of currents. The number of currents is not conserved. 



3 Analytic considerations 



The discretization of momenta via the DLCQ procedure puts severe constraints on the 
possible Fock states. It turns out that they allow for the a priori calculation of two of 
the eigenvalues. At harmonic resolution K the states 



\K) 
\K-1) 



Tr 
Tr 



{j(-i)r\ io> 

{J(-1)}*- 2 J(-2)1 |0> 



(15) 
(16) 



are unique and have K and K — 1 currents, respectively. The sectors with less than 
K — 1 currents contain more than one state, e.g. 



\K 
\K 



2)i 
2} s 



Tr 
Tr 



{J(-l)} K - 3 J(-3)} |0>, 
{J(-l)}^- 2 -V(-2){J(-l)r 2 J(-2)l |0>, 



2 < s < 



K 



(17) 
(18) 
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It is relatively straightforward to derive the expressions 

P-\K) = £{N c + N f )K\K) + OQK-l)), (19) 

P-\K-1) = l-( Nc + N f )(K-^j\K-l) + 0(\K-2)), (20) 
P~\K-2)i = + O(\K-2)), (21) 

where 0(\p)) are terms involving states with p or less currents. Let the dimension of the 
discrete Fock space be d and define fix := ^(N c +Nf)K and fi 2 '■— 2^(N C +Nf) {K — §) . 
Then the structure of the Hamiltonian matrix is 





( A 





p- = 


\ B 


C hi 



(22) 



The matrices A and B have dimensions (d — 2) x (d — 2) and 2 x (d — 2), respectively, 
and C is a real number. Clearly, fii and /i 2 are two eigenvalues of P~ . The eigenvalues 
of the mass squared operator are then 

M f = i^(l + \)K 2 and = ^(l + \)K ( K - -) . (23) 

7T 7T V 2/ 

These eigenvalues seem to diverge in the continuum limit, which would render them 
physically irrelevant. However, one can show that the eigenvalues Mf(K) lie in different 
Z 2 sectors for even and odd K and therefore cannot be connected. They rather mark 
the appearance of new states in the spectrum, as we will see. 



4 Numerical Calculations 
4.1 The Hamiltonian 

In this section we address to calculate the eigenvalues of QCD 2 by solving the eigenvalue 
problem, Eq. @, numerically. Since the operator P + is already diagonal, we have to 
construct the action of the Hamiltonian P~ on a basis stateQ Using the large iV c limit 
implies that the basis states be of the form Tr[J\ • • • Jb]|0), i.e. single-trace states. 
The current operators are subject to the Kac-Moody algebra, Eq. (|I2"D. Annihilation 
operators may thus be produced by commuting operators. The main obstacle for the 
calculations is to find a scheme to organize the terms in a convenient way. We obtain 
such a scheme by separating terms containing annihilation operators from those that 
do not. In the definition 

[A,B\ =: \A,B] + [A,B\, (24) 

\A, B] denotes the part of a commutator which consists solely of creation operators. 
Its complement \_A, B\ contains annihilation operators. It is possible to write down 

: A continuum version of the first step of this DLCQ calculation can be found in Ref. H. 
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an expression for P which involves only commutators of the first type and creation 
operators. It reads 



P-J 1 j 2 ---J b \0) 



y y j\ • • • Jki-iJki+i ■ • ■ Jk 2 -\Jk 2 +i • ■ ■ Jk P -i 

p=l k 1 <k 2 <...<k p 

x \l ■ ■ [P~, JkA , J k2 \ , ■ ■ ■ , Jk p ] Jk v+ i ■ • ■ J 6 |0>. (25) 



The commutators in this expression, and thus the Hamiltonian matrix, are linear in 
Nf by construction. To construct the Hamiltonian matrix, we have to evaluate all 
2 b — 1 commutators. In the worst case scenario we would have an exponentially growing 
number of terms (| J2p=i 2 2 ( p ~ 1 ) in the Hamiltonian, due to the |2 2 ( p_1 ) terms in a 

commutator involving p currents. Fortunately, the number of terms of leading power in 
N c grows only quadratically, like (26 — 1)6. We developed a computer code to evaluate 
Eq. (p5|) symbolically. This task exceeds typical workstation capabilities at 6 > 7. 
It is however possible to deduce the expression for arbitrary 6, because a repeated 
pattern evolves. In the large N c limit the action of P~ on a state |b;ni, . . . , n b ) := 
jjbJjl{. n i) ■ ■ ■ Jjl ( n fe)|0) with 6 currents is then 

P~\h;n x , ...,n b ) 



(26) 



+ 



+ 



+ 



~9 2 N C 
2tt 

2tt 



EE 



m-i i 



i=X \m=l ( m U i) m=l 
E _ + E — 2 H b ' n l> n2 ' • • • ' Hb ) 



b + 1; m, n 2 , . . . , rii — m, m, 



n b ) 



i=i 



11; 



m=l 



g 2 N c 6 " 1 /ni+1_1 



2tt 



E E 



1 



—1 



i=l \ m- 



~9 2 N C 



2tt 



E 



1 



— 2 |b; ni, n 2 , . . . , n; + m, n i+1 -m,...,n b ) 



m 



o( m + n ir 



"6-1 J 

^ — I |b; n 2 , ■ ■ ■ , n b -i,n b -m,n l + m) 



7,2 AT 6-2 



+ A 



2tt 



6-j 



E {(-)'£ 



i=l 



1 



n 



X |b - j; m, . . . , Tij-i, ^ n 9 , ...,n b ) 



i+j 



q=i 



v- 1 



1 



1 



. (^1 + Sg=fe_j + 1 



n b 



x 



b - j;n 2 , • • • ,n 6 _j,ni + 51 n ?) 

J E 



i=i 



i 



<J=f)-i + l ^cjj 



(E 



? =6-i+l 



n b 



6 



K 


2 


3 


4 


5 


6 7 8 9 


10 


11 12 


13 14 


states 


1 


2 


4 


6 


12 18 34 58 


106 


186 350 


630 1180 



Table 1: Number of basis states a a function of the harmonic resolution K. 



x 



3-1 b 
q=l q=b—i 

-\b-y,n j - i+2 ,n :j -. i+3 ,...,n b -i, E ra g + E n <?) 

q=l q=b—i+l 

6-2 [ b-j-1 n i+]+1 -l 



27T 



e (-y E E 



1 



1 



i=l m=0 



(m + Egti n g ) 2 (m + Eg= i+ i « 9 ) 2 



|b - j;ni,n 2 , . . . ,E n 9 + ~ m,n i+j+2 , ■ ■ ■ ,n b ) 



q=l 




E T^J |b-j;n i+ i,...,n 6 _i,^n g + n 6 ) 

Note that only the terms in lines two and five of this result contain Nf. These terms 
will be absent in the 't Hooft limit (Nf/N c — > 0) and will be dominant in the large Nf 
limit. 



4.2 Numerical results 

We are solving the eigenvalue problem, Eq. (|9]), numerically to obtain the mass spec- 
trum. Recall that we have two parameters at our disposal to study the spectra. One is 
the harmonic resolution K, which we are supposed to send to infinity. The other is the 
ratio of the numbers of flavors and colors, A = Nf/N c , of the fermions in the theory. 
The most prominent cases are A = 0, 1, oo, namely, the 't Hooft limit, adjoint fermions, 
and the large Nf limit. Note that also the unphysical parameter K might give insight 
into the spectrum, e.g. the discovery of continuum states via their characteristic K 



dependence in Ref. [ I I 



The advantages of formulating the problem with SU(N) currents rather than with 
fermionic degrees of freedom are twofold. For one, the much smaller basis allows for a 
larger resolution K and thus for more accurate results. We will also see the structure of 
the spectrum much clearer, in the sense that many, if not all, of the uninteresting multi- 
particle states are absent. Secondly, we are able to study the behavior of the spectrum 
as we couple the gauge field to different forms of matter by varying the parameter 
Nf/N c , which is an algebraic variable in the present approach. This might be used a 
tool to interpret spectra. 
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Performing the numerical calculation, we obtain exactly the same eigenvalues in 
the adjoint case as in previous works with anti-periodic boundary conditions for 



the fermions. The number of states grows exponentially with the harmonic resolution, 
cf. Table [l|. At K = 12 we are diagonalizing a Hamiltonian of dimensions^] 350 x 350, 
whereas in the fermionic approach one would have to operate on a Fock space with 4338 
states to obtain the same accuracy. To further reduce the computational effort, we can 
use the Z 2 symmetry of the Hamiltonian which is invariant under the transformation 
TJij(n) = —Jij(n). It is straightforward to convince one-self that the action of this 
operator on a state with b currents is 



T|b; m, n 2 ,..., n b ) = (-) b |Pi! ^1, Tj(n b , n b _i, . . . , n 2 )), (27) 

i=0 

where the Tj consist of Pi partial sums of the 6—1 momenta, in the sense that i runs over 
all possibilities to place 0, 1, . . . , b — 2 commas between the momenta while summing 
those momenta which are not separated by a comma, e.g. 

T\4;m,n 2 ,n 3 ,n 4 ) = |4; n 1 , n 4 , n 3 , n 2 ) + |3; n 1 , n 4 + n 3 , n 2 ) 

+ |3; n 1} n 4 , n 3 + n 2 ) + |2; n h n 4 + n 3 + n 2 ). (28) 

In order not to complicate the construction of the Hamiltonian, we determine the Z 2 
parity of an eigenstate a posteriori by calculating the expectation value of the operator 
T in this state. The separation of the Z 2 odd and even eigenfunctions is useful when 
interpreting the results, because it reduces the density of eigenvalues to roughly a half. 



4.3 The 't Hooft limit 

To test the consistency of our approach, we compare to the well-known results in the 
large N c and the large Nf limit. We find complete agreement. First let us consider 
the large N c (or 't Hooft) limit, where we should recover the results of the 't Hooft 
model [0]. The spectrum of large N c QCD in two dimensions with massless fundamental 
fermions asymptotically has the form M 2 = 7i 2 n for large integer n, where the mass 
squared is in units g 2 N c /ir. We find that at K = 14 the lowest ten single-particle states 
have masses^ 

M 2 = 5.88, 14.11, 23.04, 32.27, 41.68, 51.24, 60.93, 70.76, 80.97, 90.90, (29) 

which is in very good agreement with 't Hooft's numerical solution [|TJ. The actual 
spectrum, Fig. is a mixture of single- and multi-particle states. The multi-particle 

2 It is not clear a priori how to construct a Fock basis for the current operators. A naive expectation 
is that single-trace states modulo cyclic permutations provide such a basis. We convinced ourselves that 
this is actually the correct choice by explicitly performing a large N c Gram-Schmidt orthonormalization 
for small K on these states. 

3 We performed the continuum limit here and in the sequel by fitting the data to a polynomial of 
second degree in 1/K. 
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Figure 2: Spectra of the 't Hooft limit in the Z 2 even (a) and odd (b) sectors. Solid lines 
connect associated eigenvalues as K varies. Short dashed lines connect analytically 
calculable eigenvalues. Long dashed lines are extrapolations towards the continuum 
limit. Masses are in units g 2 N c /n. 



states decouple completely from the single-particle states. They appear here, because 
we do not use an orthonormal basis. We have performed a calculation with an or- 
thonormal set of states up to K = 6 and found that then only single-particle states are 
present. We were able to identify all multi-particle states as composites of two or more 
single-particle states. The masses M 2 of multi-particle states are given by 

M 2 mp (K) _ Ml(n) M^(K-n) 

~^T~ ~ + K-n ' m 

where Mp.(n) are the masses of single-particle states at harmonic resolution n. These 
masses are exactly reproduced in the spectrum. This identification in turn made it 
possible to detect the single-particle states hidden amongst the multi-particle states. 
Surprisingly, the eigenfunctions have no significant structure, apart from the special 
form of the lowest state in each Z 2 sector. In particular, we were unable to distinguish 
single- from multi-particle wavefunctions by their shapes. 



4.4 Adjoint fermions 

Consider now the adjoint spectrum, Fig. |3|. If we look at the eigenvalue trajectories 
(mass squared as a function of K), the structure of the spectrum looks similar to the 
't Hooft case. We see immediately three single-particle candidates which qualify by 
their straight, smooth trajectories. In the continuum limit they have the eigenvalues 

M B = 10.84, M| = 25.73, Ml = 45.66. (31) 
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0.1 0.2 0.3 0.1 0.2 0.3 

1/K l/K 



Figure 3: Spectra of the theory with adjoint fermions in the Z2 even (a) and odd (b) 
sectors. Solid lines connect associated eigenvalues as K varies. Short dashed lines 
connect analytically calculable eigenvalues. Dashed lines are extrapolations towards 
the continuum limit. Masses are in units g 2 N c /ir. The dash-dotted trajectory is the 
analog of the third single-particle state in the 't Hooft limit. 

The expectation for the asymptotic solution is here |J 

Ml, m ,..., nk = 27T 2 (m + n 2 + • • • + n k )- n G 2Z, (32) 

again in units g 2 N c /ii. The single-particle masses appear to be at roughly twice the 
values of the 't Hooft limit, whereas the low-lying 'multi-particle f] masses stay at more 
or less their 't Hooft values. The crucial difference between the spectra is that the 
'multi-particle' eigenvalues of the adjoint case coincide only approximately with the 
values from the mass formula, Eq. fl3U|). They are systematically lower than expected. 
As such, these states appear to be loosely bound composites, i.e. they have to be 
interpreted as single-particle states, at least at finite K. We can see that these states 
are a vital part of the spectrum by comparing the eigenvalues at M 2 ~ 32, K = 13 
in Fig. |(a) with those of Fig. |(a) at M 2 ~ 42, K = 12. In the latter case the 
eigenvalues repel each other, leading to a deformation of the eigenvalue trajectory of 
the (conjectured) single-particle state B 3 . This means that the 'multi-particle' states 
are interacting! By varying the parameter A, we were able to follow this deformation. 
As Nf grows, the lowest multi-particle states of the 't Hooft spectrum 'move through' 
the single-particle state with mass squared M 2 = 23.04 in the 't Hooft limit, and 
produce these deformations in a trajectory. In the 't Hooft limit, Fig. @(a), we see 

4 We use quotation marks here and in the following because it is not a priori clear that these states 
are indeed multi-particle states. 
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Figure 4: Adjoint fermions: Left (a): The current number distribution functions of 
selected states at different harmonic resolution K. The lower three rows are distribu- 
tions of (conjectured) single-particle states. The state in the top row is not present at 
K — 11. Right (b): Wavefunctions of selected eigenstates at K — 14 as functions of 
the basis state number divided by the total number of states (1180). The dashed lines 
mark points of changing parton content (number of currents in parenthesis). 



no distortion of the single-particle trajectory though the corresponding eigenvalue is 
almost degenerate with a (decoupled) multi-particle stateQ. 

The question arises, whether these loosely bound states become multi-particle states 
and decouple in the continuum limit in the adjoint spectrum. For the states between 
M 2 = 20 and M 2 = 40 in the Z2 even sector, the answer is certainly yes. We performed 
a fit to the data finding that the eigenvalues converge towards a single point, namely 
the two particle continuum threshold at M 2 = 22.85 = 4Mf. , i.e. twice the mass of the 



lowest fermionic state |l.TJ. The deviations from the mass formula, Eq. ([Hj), vanishes 
faster than -^p, where (3 > 2. 

The deformation of a trajectory renders it questionable if the associated state is 
really a single-particle state. As K grows there will be more and more 'multi-particle' 
states to distort the trajectory. On the other hand, the 'multi-particle' sector will in- 
teract less and less as K increases. The number of (low-lying) 'multi-particle' states 
grows linearly with K. If we can take the deviation of the 'multi-particle' states from 
the mass formula (proportional to jrp) as a measure of the coupling, we have to con- 
clude that the single-particle trajectory will be less and less distorted and thus a true 
single-particle state will emerge as K — > 00. Although the first two points of the tra- 



5 The same deformation occurs in the fermionic spectrum, namely in the trajectory of the second 
Z2 even single-particle state of Ref. jfl| . 
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Figure 5: Left (a): Spectrum of the Z 2 even sector in large Nf limit. Note that masses 
are in units g 2 Nf/ir. The dashed line is the extrapolation of the lowest eigenvalue 
towards the continuum. Right (b): Mass squared eigenvalues as functions of A = 
Nf/N c . Masses are in units g 2 N c /7r. 

jectory at K = 4 and K — 5 are multi-particle states with fermionic constituents^, 
we feel quite safe about the single particle nature of this state, because we find that 
the same situation occurs in the case of the lowest Z 2 odd trajectory, which represents 
undoubtedly a single-particle state, due to its complete isolation in the spectrum. To 
check whether all 't Hooft mesons become adjoint single-particle states, we followed the 
next two Z 2 even single-particles of the 't Hooft spectrum as we varied A in small steps. 
One is shown as the dashed-dotted line in Fig. |3] just above the state converging to 
M 2 =45.99. Their irregular shape makes it unlikely that they are single-particle states. 

The eigenfunctions and the probability distributions in parton number were helpful 
to identify states at different K in the fermionic formulation of the theory |Q, Here, 
with a non-conserved parton number, these quantities will be of limited usefulness. 
In the fermionic formulation the single-particle states have a high probability to be 
states of a definite parton number [[J. In the present work, we find this behavior only 
for the highest states in each Z 2 sector, which are states with K or K — 1 currents, 
cf. Sec. ^. Note, however, that the conjectured lowest single particle states B\, B 2 , and 
B3 have a large two-current contribution, no two-current, and no two- and three-current 
contribution, respectively. We definitely cannot use information about the probability 
distributions to distinguish states: if we have a close encounter of two eigenvalues in 
the spectrum, one of which is presumably a single- and the other one a multi-particle 
state, we find that the probability distributions are almost identical. For the higher 
states significant structures in the parton probability distributions are missing alto- 



6 One state consists of Mf. (.fif=§)=5 and M\ (-K"=|)=12.5, the other of the first mass plus 
(K=l)=9. 
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Table 2: Eigenvalues in the adjoint case, the 't Hooft limit and the large Nf limit. 
B\ and B 2 are the lowest lying single particle states in the Z 2 even and odd sectors, 
respectively. The eigenvalue in the large Nf limit is actually the continuum threshold. 
Note that the masses are given in units g 2 N c /ir in the first two cases and in units 
g 2 Nf/n in the third. 

gether. Concerning the wavefunctions, we find that the lowest eigenstate has a peak in 
the two-particle sector, cf. Fig. |](b), and has some half-sine structure in higher current 
sectors. We find no obvious structure in other states that could give a hint of how 
to construct solutions of the theory analytically. However, there might be hints of an 
additional structure in the eigenfunctions: we found several wavefunctions which have 
a noticeable amplitude drop at the boundaries of parton number sectors. As an exam- 
ple, consider the drop of the wavefunctions of the 'single-particle' states (B 1 , B 2 , B 3 ) at 
the boundaries of the 2, 3, and 4 current sectors, respectively, in Fig. £|(b). In general, 
'single-particle' states exhibit a sharp change of the sign of the amplitude at the bound- 
aries of current sectors. The 'multi-particle' state (M 2 = 22.03) does not exhibit such a 
behavior. In principle, the wavefunction should factorize if it is a continuum state con- 
sisting of two non-interacting bound states. Disentangling wavefunctions along these 
lines has been attempted in Ref. [BTJ; here, however, this task seems intractable, due 



to the non-conserved parton number. 



4.5 Large Nf limit and intermediate cases 

When going over to the large Nf limit, A — > 00, the operator P~, Eq. fl2"6]) , has a dom- 
inant kinetic term. In the two-particle truncation the spectrum is purely kinematical. 
In the full theory a residual interaction is present. In the spectrum, Fig. |5|(a), we see 
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a continuum of states starting at M 2 = 4g 2 Nf/ir. The only particle in the spectrum 
is thus a meson with mass M 2 4 = g 2 Nf/ir. This is in full agreement with results from 
analytic calculations [|7|]. The Z2 odd sector looks similar, with the continuum starting 
at M 2 = 9 = 3 2 M|, f . The single-particle state is absent in our calculations due to the 
tracelessness of the currents, i. e. the absence of the one current state. 

Let us look at the eigenvalues as a function of Nf/N c . Fig. j|(b) shows the first[] 
calculation of the spectrum of two-dimensional Yang-Mills theories as a function of 
the continuous parameter Nf/N c . We see that the (low-lying) trajectories are linear, 
as is expected from the analytic calculations in Sec. || The parameter Nf actually 
plays a role similar to that of a mass term in the 't Hooft model |18|]. Indeed, the 
curves obtained by taking the continuum limit of the single-particle trajectories in 
Fig. [5](b) look very similar to results of massive two-dimensional QCD |[20|| . Other 
DLCQ calculations of two-dimensional massive QCD ^] were performed at finite 
iV c . In particular, it was shown in Ref. [^TJ that the spectrum of massive QCD in two 
dimensions is well described by the large N c approximation. Deviations from the large 
N c behavior, found typically at small fermion masses m where the relevant coupling 
g 2 N c /m 2 is strong, might just be artifacts of the DLCQ approach. This is in line with 
results of Ref. [|l^], showing that also in adjoint QCD in two dimensions the 1/N C 
corrections are very small; a fact that still awaits proper explanation. 

We see a lot of level crossings in Fig. |(b). The fact that the trajectories do not 
intersect is a finite K effect. We note that the plot shows a smooth behavior at A = 
1. There is no hint that the adjoint case could be a special point in the parameter 
space. The eigenfunctions in the large Nf limit look very different than those found 
in the adjoint and the 't Hooft case. Typically, the amplitude vanishes, except for a 
few delta-like peaks, reflecting the fact that the spectrum of this theory consists of 
continuum states of non-interacting mesons. For comparison with other work, we list 
the eigenvalues obtained in the numerical calculations in the 't Hooft limit, with adjoint 
fermions and in the large Nf limit in Table |2|. 

5 Summary and discussion 

In this note we studied two-dimensional QCD using SU(N C ) currents as basic degrees 
of freedom. This enabled us to calculate the Hamiltonian matrix and the spectrum of 
any two-dimensional SU(N C ) gauge theory coupled to an arbitrary number of fermions. 
Working in the DLCQ framework, all relevant quantities become functions of the pa- 
rameter A = Nf/N c and of the harmonic resolution K. We constructed the Hamiltonian 
matrix explicitly in the Fock basis of single-trace states. 

We reproduced all known results in the A parameter space: in the 't Hooft limit 
(A — > 0) we obtained the well-known linear spectrum, in the adjoint case we find exactly 
the eigenvalues of previous works 0, [13], |TTJ, and we identified the single meson of the 
large Nf limit. Moreover, in the case of adjoint fermions we were able to confirm the 



single-particle nature of the lowest Z2 odd boson, which was hidden in Ref. |TT[ among 



7 M. Engelhardt |17[ presented a truncated version of this calculation. 
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(trivial) multi-particle states which are absent in our calculation. We provided evidence 
for the single-particle nature of the boson B% at M 2 =45.99, which makes the existence 
of other single-particle states above the continuum threshold at M 2 = 4MJ. very 
likely. We were able to explain the deformations in the trajectories of the eigenvalues 
as repulsions of eigenvalues. Unfortunately, the amount of remaining 'multi-particle' 
states in the spectrum makes it impossible to decide whether or not the identified single- 
particle states form an infinite Regge trajectory, or if even a multi-Regge structure 
exists, as was suggested in Ref. [fT9| . What we were able to do is to eliminate all states 
from the list of single-particle candidates which do not appear in our calculations, but 
in the bosonic sectors of previous works. By construction, our approach contains all 
single-particle states ||. 

In Ref. [II] it was suggested the massless limit is reached only in the continuum 
limit K —>■ oo. We work in a manifestly massless approach and obtain exactly the 
same eigenvalues for all K. It would be very interesting to repeat the calculation of 
Ref. [II] for finite m, as was done with different motivation in Ref. |I3 |. With our 



results at hand, one can focus on a much smaller set of single-particle candidates, and 
see how these states develop as the mass is turned on. Our guess is that this transition 
is continuous. An interesting related question is how to distinguish screening from 
confinement when only the mass spectrum is known. 

We operate at higher numerical precision than previous work. This allowed us to 



prove numerically that the continuum threshold found in Ref. [I !j is indeed exactly at 
four times the mass squared of the lightest fermion. We note however that the inter- 
pretation of the continuum states is not completely clear. We find that the deviation of 
their masses from the expected free many-body masses, Eq. fl30|), vanishes faster than 
1/K 2 as we go towards the continuum. If these states form a continuum at exactly 
the expected threshold, one has to conclude that they decouple completely from the 
single-particle states and that their coupling is an artifact of the finite resolution K. 

To summarize, we presented a refined and quite coherent picture of two-dimensional 
QCD by using a new computational tool. We hope that this approach will prove 
powerful and that it will yield new qualitative insight, too. We mainly focused on 
quantitative improvements in the present work. We hope that this work will provide 
sufficient input for future enterprises to understand this theory, which shares some of 
the key features of full QCD, much better. 
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